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Abstract 

For Klein-Gordon equation a consistent physical interpretation of wave 
functions is reviewed as based on a proper modification of the scalar product 
in Hilbert space. Bound states are then studied in a deep-square- well model 
where spectrum is roughly equidistant and where a fine-tuning of the levels is 
mediated by VT — symmetric interactions (composed of imaginary delta func- 
tions) which mimic creation/annihilation processes. 
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1 Klein-Gordon equation 



1.1 Pseudo-Hermitian Feshbach-Villars Hamiltonian 

As long as the most common relativistic Klein-Gordon (KG) operators are partial 
differential operators of the second order with respect to time, the time evolution 
of the wave functions ty( KG >(x, t) must be studied together with their first time 
derivatives i dt ty( KG ^ (x, t). After the routine Fourier transformation we arrive at the 
Feshbach-Villars (FV, [I]) non-Hermitian eigenvalue problem 



m = E\ip) , hw = u I (i) 




where the two wave-function components may be marked as \D} (="down compo- 
nent") and \U} (="up component"). For the description of the bound states in one 
dimension the two-by-two partitioning in allows us to extract \U} = E \D} and 
to replace our Klein-Gordon equation by its reduced form 

U KG) \D n } = e n \D n }, n=l,2,... (2) 

with squared energy E 2 abbreviated as e and with the "large" Hilbert space Ti of kets 
\ip) reduced to the "smaller" Hilbert space TC( C ) of the curly-ket "down" components 

\D n } m 



1.2 Biorthogonal bases 

The "right" eigenkets \D n } will not carry all information about h( KG ^ whenever 
^(-KG)] ^ \h( KG )^ . Then, the parallel Schrodinger-type problem generates different 
eigenkets marked by the double curly ket symbol. The latter sequence may be re- 
read as the left eigenvectors of our original operator h^ KG \ related to the same (by 
assumption, real) eigenvalues e n = 

{{L n \U KG ^ = K 2 n {{L n l n = l,2,.... (3) 
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It is well known that the set of the bras {{£„ | and kets \D n } is bi-orthogonal [2], 

{{i m |A,} = (° Jn for 

{ Qn T { m = n, 

and that it forms, usually, a basis in the infinite-dimensional Hilbert space TCm. 
Then, we may decompose the unit operator and/or derive the bi-orthogonal spectral 
representation of the Hamiltonian in T~lr c ), 

oo i oo 2 

Ac) = E I A.} ~ {{Ln\ , U KG > = £ \D n } {{L n \ . (4) 

n=l ^ n=l 6n 

The overlaps Q n need not be all of the same sign. 



2 Relativistic observables 
2.1 G— quasi-Hermiticity 

In the space H = H( c ) © W( c ) °f the eigenstates of H^ FV ^ we have to consider the pair 
of conjugate equations 

H{FV) | n ( ±)) = |n ( ±)); ((n (±)| ^(jto = ±Kn ((n (±)| (5) 
Both the left and right eigenstates have the two-component structure, 
| m (±))>=( ), | n (±)) = f ±K » |jD "> 

y ±/« m |L m }} y y |D n } 

and form the bi-orthogonal set in the "bigger" space 7i, 

((m {v) \rS v,) ) = 5 mn 5 uu/ ■ /i^ , /i^ = ±2/t m £ m , i/, z// = ±1 . 
It is expected to be complete and useful, 

OO 1 

1= E E l" (r> >-^«n (T) l, (6) 

t=±1 n=l f^n 
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hw= e EI« Ct, >^«« (t) I = E 

t=±1 n=l /in n=l Z " n 

Let us now assume that at a given H^ FV \ equation 

ff^f^ff^rT 1 (7) 

possesses a positive and Hermitian solution r] + = 9 > 0. Such an operator may play 
the role of a metric and induces the following specific scalar product in 7i, 

(|V>i> © l^ 2 » = © l^ 2 > = (^1^2)0,^^), IV>i> e W, 1^2) e W. (8) 



(|n«) ((n«|) + (|n^> ((n^| 



This product generates the norm, = \J(ip\ip) (physical)- In terms of the later 

product and metric we may call all the operators A with the property = A @ _1 
quasi- Hermitian and treat them as observables (see [3] for a deeper outline of some 
more sophisticated mathematical details). Indeed, we have 

(|^l)©|AV>2>) = (|AV>l)0|^2>) (9) 

so that the probabilistic expectation values A (physical) are mathematically un- 
ambigously defined. 

2.2 Explicit constructions of the metric O 

Let us assume non-Hermiticity of the type h^ KG ^ 7^ [h^ KG ' 1 ]^ = Vb^ KG ^V in the 
smaller space Hi a (here, V is operator of parity). Then, a consistent physical meaning 
may still be assigned to all the relativistic bound states, provided only that in the 
bigger space TC we find a suitable physical metric 0. For this purpose we may employ 
the ansatz 

00 

©= E E \n {T) ))M^((m^\, 

t.t/=±1 m,n=l 

the backward insertion of which in (JJJ) gives the condition 

TK n M^=M^T!K m 
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with the set of solutions M^p = u>^8 nm 5( TT/ ) numbered by the free parameters 
uj^\ The Hermiticity and positivity constraints restrict the freedom of the choice 
of both the optional sequences cu^ to the real and positive values, uj^' > 0. Vice 
versa, any choice of the latter two sequences defines an eligible operator of the metric 

oo 

© = 6*00= E E |n«»o;W((nW|. (10) 

T=±ln=l 

Its inverse 

OO 1 

r=±l n=l UJn \ [In \ 

is similar. In terms of the metric 0, the formal bound-state wave functions re- acquire 
the standard probabilistic interpretation. 

3 Models with complex point interactions 

In a way inspired by the success of several non-relativistic studies of VT— symmetric 
models with point interactions [3] and by the encouraging experience we made in our 
paper [5] we shall combine the infinitely deep square-well real part of the potential 
[V(x) = oo for all x ^ (—1,1)] with the following purely imaginary delta-function 
formula for its remaining part, 

c 

V( x ) = -i&5(x + a t )] , xe(-l,l), (12) 

at real couplings & and ordered points < a± < a 2 < . . . < ac-i < ac < 1. 

3.1 Wave functions 

The key advantage of our V(x) in (|12|) is that the VT— symmetrically normalized 
coordinate representants ip{x) = ip*(—x) of \D} in eq. (J2J) remains piecewise trigono- 
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metric. At each real and positive bound-state energy e = k, 2 we shall have 

^ , L ) ( a; ) = ( a £ -i>Pt) sin K(a e+1 + x) + (7^ - z<^) cosK(a m + x), 

x G (— a^+i, -at), 



ip(x) 



X) 



/j, cos kx + i v sin /cx, 



x G (— ai, ai 



= ( a £ + * A) srn K ( a e+i — x ) + (t« + i $l) cos K ( a m _ x ), 

x G (a^,o <+ i), 

ipR \x) = (ac + i (3c) sin k(1 — x), x E (ac, 1), 1 < £ < £. 
Its differentiation as well as continuity conditions 

^t X \~o-i) = *f{-a.i), l = £,£-l,...,2, 

4 0) (-«i) = ^(01) = 

V^ +1) (o/+i) = £ = 1, 2, 1, 

enter the definition of the action of the delta functions, 

e = cx-i 



-at, 



4°> 1 



-Oi) 



-o x . 



-a £ j, £ 



,2, 



-01. 



4 0) («i; 



1. 



(13) 



(14) 



(15) 



^=1,2,...,£ 

After the insertion of the ansatz (|13j) . the set of formulae (|14j) and (J15)) may be read 
as a homogeneous linear algebraic system of 4£ equations for the 4£ unknown wave- 
function coefficients ac, (3c, ■ ■ ■ ,v- The secular determinant T>(k) of this system 
must vanish so that the not too complicated transcendental equation 



V(k) = 

determines finally the set of the bound-state roots k = K n at n 



(16) 



1,2, 



3.2 Energies at the simplest choice of C = 1 

At £ = 1, potential (JT2~j) degenerates to the most elementary double- well model with 
the single coupling £1 = £ and one displacement a% = a jHj- Out of the related eight 
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real constraints (|14|) and (|15|) only four are independent and define the four real 
coefficients a,\ = a, (3\ — (3 and /x and v as an eigenvector of a four-by-four matrix 
with the secular determinant 




Numerically, the first term would give us the well-known square-well spectrum at 
£ = 0, the completeness of which is controlled by the Sturm-Liouville oscillation 
theory jjj. As long as all the roots n n = K n (£) are smooth and real functions of £ at 
the smallest couplings, n n {0 ~ wr/2 + 0(£ 2 /n), our explicit construction confirms 
the general mathematical prediction jS] that the influence of the non-Hermiticity will 
be most pronounced at the lowest part of the spectrum. 



3.3 The next choice of C = 2 

In the quadruple-well potential (|12|) with £ = 2 we may shorten b and 

drop the two redundant subscripts in 71 = 7, Si = 5. In the eight- dimensional matrix 
of the system the elimination of four unknowns is either trivial [7 = a 2 sin«(l — b), 
5 = /3 2 sin/«(l — b)] or easy [ai = cti(ct2, P2), Pi = Pi{&27 fa)]- We end up with a 
four-by-four matrix problem and with the secular determinant 

V{k) = Z> (0) (k) +£>&)(«) +2? (6) («) + (18) 

1 £ 2 
^(o)(«) = - » sin 2k , = sin 2Ka i " sin 2 [fi;(l - a,-)], j = 1, 2, 

^(5i6)( K ) = — I ~^T~ SU1 + -^j^ sin 2 [fi;(6 — a)] 1 sin 2 [/s:(l — b)] . 

This secular determinant correctly degenerates to the previous C = 1 formula in 
both the independent limits of £1 — > and £ 2 — > 0. 
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3.4 Simplifications at the rational aj 

Let us return to the secular eq. (|T7|) with C = 1 and choose a = 1/2 [6 . This leads 
to a factorization of T>(k) and to the pair of the eigenvalue conditions 

cosK 2m _! = — — , sm/c 2m = 0, m = l,2, ... (19) 

with the second series of equations being exactly solvable, K 2m = m7r. 

At the next choice of a = 1/3 we factorize eq. ()17p in the similar manner and get 
the series of the £— dependent roots specified by the implicit definitions 

cos 3** = ^2 _ 4Ag 2 » P = 1> 2, 4, 5, 7, 8, 10, . . . (20) 

complemented by the closed formula for all the skipped roots of the factor sin2«;/3 
which remain £— independent and read K% m = 3rrm/2 with m = 1, 2, . . .. The regular- 
ity of such a pattern of the £— independent roots is easily prolonged to the decreasing 
sequence of a with ft 4m = Inm at a = 1/4 and all m — 1, 2, . . ., etc. 

The less elementary composite choice of a = 2/3 may be observed to give the 
same factor as at a = 1/3 and, hence, the same £— independent series of the roots 
^3m = 3m7r/2 with m = 1,2,.... The implicit formula for the remaining roots 
is a slightly more complicated quadratic equation in the trigonometric unknown 
X = cos2k/3, 

(4k 2 -£ 2 ) X 2 + £ 2 X-k 2 = 0. (21) 

Its trigonometric part X may be eliminated in the form resembling eq. ()20j). 

One of the important consequences of the existence of the elementary formulae 
for the rational a is that they allow us to perform an elementary analysis of the 
qualitative features of the n— th root K n during the growth of the strength £ of the 
non-Hermiticity. During such an analysis one discovers that these levels are either 
"robust" (marked by a superscript, and remaining real for all £) or "fragile" 

(such a will merge with another nffl at a "critical" £^ while the pair will 



S 



complexify beyind this "exceptional" J2| point). For illustration let us display this 
pattern in the three simplest spectra, 

JF) OR) (F) (R) (F) (R) n -Ml 

JF) (F) (R) (F) (F) (R) n-M-X 
ft 1 , ft 2 , ft 3 , ft 4 , ft 5 , ft 6 , . . . , a — i/o 

K GF) „(*■) JR) JF) (F) (R) (R) (R) (F) _ 

ft x , ft 2 , ft 3 , ft 4 , ft 5 , ftg , ft 7 , ftg , ftg , ft 10 , ft-11 , • • • , U — 
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